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ABSTRACT 

We present the results of a broadband spectral study of the transient Low Mass 
X-ray Binary Aql X-1 observed by Suzaku and Rossi X-ray Timing Explorer satellites. 
The source was observed during its 2007 outburst in the High/Soft (Banana) state 
and in the Low/Hard (Extreme Island) state. Both the Banana state and the Extreme 
Island state spectra are best described by a two component model consisting of a 
soft multi-colour blackbody emission likely originating from the accretion disk and 
a harder Comptonizcd emission from the boundary layer. Evidence for a hard tail 
(extending to ^ 50 keV) is found during the Banana state; this further (transient) 
component, accounting for atleast ~ 1.5% of the source luminosity, is modeled by a 
power-law. Aql X-1 is the second Atoll source after GX 13-1-1 to show a high energy 
tail. The presence of a weak but broad Fe line provides further support for a standard 
accretion disk extending nearly to the neutron star surface. The input photons for 
the Comptonizing boundary layer could either be the disk photons or the hidden 
surface of the star or both. The luminosity of the boundary layer is similar to the disk 
luminosity in the banana state and is about six times larger in the extreme island 
state. The temperature of the Comptonizing boundary layer changes from ^ 2 keV in 
the banana state to ^ 20 keV in the extreme island state. 
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1 INTRODUCTION 

Low mass X-ray binaries (LMXB) containing weakly mag- 
] netised neutron stars (NS) are classified into Z and Ato ll 
$^ sources (|Hasinger &: van der Klislll989l :l van der Klisll2006h . 
' The classification is based upon the different tracks these 
two types of sources trace in the colour-colour diagrams 
(CD) and their associated timing properties in the different 
spectral states. Even though detailed timing and spectral 
differences between the two class of sources are not com- 
pletely explained, higher accretion rate is generally believed 
to characterise the Z as compared to the Atoll sources. Atoll 
sources have lower luminosities ( 0.001-0.5 LEdd) than Z 
sources and they have fragmented CD consisting of banana 
state (BS) and Extreme island state (EIS). There is also a 
transitional state between these two states referred to as Is- 
land state. In the BS the source luminosity is high and the 
observed energy spectrum is soft (high/soft state). On the 
other hand in the EIS the source luminosity is low and the 
observed energy spectrum is hard (low/hard state). 

Irrespective of the type of source, there are two dif- 
ferent approaches to model their energy spectra, both con- 
sisting of a soft/thermal and a hard/Comptonized compo- 
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nent f see iBarretl 1200 ll for a review). In one approach the 
accretion disk is considered to be emitting the soft ther- 
mal component while the boundary layer emits the hard 
Comptonised component ("Eastern Model" iMitsuda et al.l 
[l989. i In the other approach, the boundary layer emits the 
soft thermal component and the inner hot accretion disk 
emits the hard C omptonised component of the spectrum 
("Western Model" IWhite. Stella, fc Parmailll988l V To dis- 
tinguish between these two approaches good quality broad- 
band spectral data is required. Beppo-SAX broadband (0.1- 

200 keV) observations show that Eastern-like models are 
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pre ferred for several NS-LMXB (e.g. [G uainazzi et al. 199g; 

Lav agetto et al.|[20o3 : iTarana. Bazzano & Ubcrtini 2008). 

Suzaku observations cover a broad energy range of 0.7- 
70 keV with good spectral resolution. This capability of 
Suzaku allows to distinguish between several modelling ap- 
proaches for NS-LMXBs 
the Z-source LMC X-2 (jAgrawal 

like, as in the case of the pe riodic burster XB1323-619 



Eastcrn-likej_ as in the case of 



Mlsrall2009l '): Western- 



iBalucihska-Church et al"]|2009l ): scattered- in dust emission 
plus direct NS and /or disk blackb ody emission, as fo r the 
eclipsing LMXB SAXJ1745. 6-2901 (|Hvodo et al.ll2009l ). 

Aql X-1 is a recurrent soft transient LMXB with 
an o ptical companion sta r of V = 21.6 K7V spectral 
type (|Chevalier et al.lll999l ) and orbital period 18.9 hours 
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Table 1. Suzaku observation details 



Suzaku 


Suzaku Obs 


Suzaku Obs 


Source count rate 


Corresponding RXTE 


RXTE Obs 


Obs Id 


Start time (MJD) 


Exposure (s) 


from Suzaku obs(cnt/s)t 


Obs Id 


Start time (MJD) 








XIS 


PIN 






402053010 (BS) 


54371.64 


13825 


139.90 ±0.18 


1.13 ±0.01 


None 




402053020 (Ell) 


54376.99 


15132 


13.27 ±0.02 


0.79 ±0.01 


93405-01-04-01 


54376.98 


402053030 {EI2) 


54382.21 


19711 


15.77 ±0.02 


0.96 ±0.01 


93405-01-05-01 


54382.34 


402053040 {EI3) 


54388.33 


17915 


12.10 ±0.02 


0.79 ±0.01 


93405-01-06-00 


54388.36 



t Errors quoted on count rates are 1 sigma errors 



IWelsh. Robinson, fc Youngj |2000| ). The distance to the 
sourc e is estimated to be anywhere between 4 kpc to 6.5 
kpc l|Rutledge et al.l l200ll ) . The recurrent transient X-ray 
outbursts observed are thought to occur due to the rmal in- 
stabilities in the accretion disk (Ivan Paradiis|[l996l l. In Aql 
X-1 the typical outburst recurre nce time is 200 days, with 
average durations of 40-60 days jSimonl [2002i V The source 
also shows Type-I X-ray bursts which are believed to be run 
away thermonuclear b urning of the matter accreted onto 
the surface of the NS (|Lewin. van Paradiis fc Taamlll995l : 
iBildster] ll998l V High frequency quasi periodic oscillations 
(QPOs) between 750-830 Hz a nd burst oscillation s have also 
been discovered in this source (jZhang et al.lll99^ ). Recently 
an episode of coherent pulsations at 550.27 Hz lasting for 
150s was discovered during the peak of the 1998 outburst 
ijCasella et al.|[200^ ). The occurrence of Type-I bursts, burst 
oscillations and the coherent pulsations in Aql X-1 confirm 
that the compact object is a rapidly spin ning neutron star. 
iReig. van Straaten. fc van der Klia (|2004l) analysed five X- 
ray outbursts of Aql X-1 observed using RXTE-PCA be- 
tween 1997 February to 2002 May. They conclude from their 
timing and spectral studies that this is an Atoll source. 

Here we study the broad band spectrum of Aql X-1 
using the pointed Suzaku and Rossi X-ray Timing Explorer 
(RXTE) observations. 



2 OBSERVATION AND ANALYSIS 

Aql X-1 went into outb urst during August 2 007 and was 
observed using Suzaku (|Mitsuda et all |2007| ) and RXTE 
l|Jahoda et al.lll996l ). The source flux increased from ~ 10 
mCrab to ~ 0.2 Crab during the outburst before it finally 
decays to a flux of ~ 5 mCrab. Figure [T] s hows the one 
day a veraged RXTE All Sky Monitor (ASM; iLevine et al] 
1 19961 ) light curve of this outburst. Seven pointed observa- 
tions are m ade using Suzaku. Th e X-ray Imaging Spectrom- 
eter (XIS; iKovama et al] |2007| ) was operated in the 1/4 
window mode during all observations. An additional 0.5s 
burst option was used during the flrst observation (Obs 
Id 402053010). Since all these observations were done after 
November 9, 200 6 no XIS2 data are av ailable. Hard X-ray 
Detector (HXD: pTakahashi et all |2007^ observations were 
done with XIS nominal pointing. For the last two Suzaku 
observations the count rates are low and the spectra ob- 
tained from XIS and HXD have poor statistics. The flfth 
Suzaku observation data could not be analysed since the 
attitude file provided by the HEASARC site is corrupted. 
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Figure 1. Figure shows the outburst light curve of Aql X-1 
as observed using the RXTE-ASM (1.5-12 keV). The square 
points marked at the top of the figure show the times of 
RXTE-PCA observations and the circled points below the 
squares show the times of Suzaku observations. We note here 
that a flux of 1 C rab corresponds to ~ 75 ASM counts/s 
l|Levine et al.lll996h . 



Hence these three observations are not used for the work 
presented here. For more details see Table [T] 

There are 16 pointed RXTE Proportional Counter Ar- 
ray (PCA) observations between 54360 MJD and 54400 
MJD, seven during the peak of the outburst and nine after 
the decay of the outburst (Figure [TJ. Three of the PCA ob- 
servations are simultaneous to the second, third and fourth 
Suzaku observations (Table [T]). In the following subsections 
we give the analysis details of data reduction for RXTE PCA 
and Suzaku respectively. 



2.1 RXTE-PCA 

PCA data are used to create the CD of Aql X-1 for the 2007 
outburst. To normalise the Aql X-1 colours by Crab colours 
we used the three PCA Crab observations between 54360 
MJD and 54400 MJD. Background subtracted lightcurves 
of Aql X-1 and Crab were extracted using Std2 mode data 
in energy bands 3.0-4.5 keV, 4.5-6.0 keV, 6.0-9.7 keV, 9.7- 
16.0 keV with 16s time resolution. To eliminate artifacts 
that may arise due to different gain corrections for the dif- 
ferent PCUs only photons collected in the third PCU were 
used since this PCU was on during all the observations. The 
Soft colour (SC) is defined as the ratio of count rates in en- 
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ergy bands 4.5-6.0 keV and 3.0-4.5 keV and the Hard colour 
(HC) as the ratio of count rates in energy bands 9.7-16.0 keV 
and 6.0-9.7 keV. Average Crab colours obtained are SC = 
1.986 ±0.002 and HC = 0.5992 ±0.0007. FigureHis the CD 
obtained by plotting normalised Aql X-1 HC versus SC. All 
normalised Aql X-1 colours having error greater than 5% of 
the observed colour were ignored. This amounted in rejecting 
14 out of 307 points of the CD. From Figure [2] we conclude 
that during the second, third and fourth Suzaku observa- 
tions, Aql X-1 was in the EIS (compare Figure [2] of this 
pape r with Figure 2 of iReig. van Straaten. fc van der Khd 
I2OO4I ) We will refer to these three Suzaku observations as 
Ell, EI2, EI3 respectively. It is noted that EI2 observation 
has the best statistics in the EIS and also EI2 is brighter 
than Ell indicating irregular decline in the lightcurve (Ta- 
ble [T]). 

There are no simultaneous PCA observations during the 
first Suzaku observation to ascertain the spectral state of 
the source. Hence we use a different method to find the 
normalised Aql X-1 SC and HC during this observation. 
The best fit spectral model that can reproduce the observed 
Suzaku spectrum is used to simulate a PCA spectrum. The 
f akeit command of XSPEC package is used for this simula- 
tion. Average count rates in the four energy bands required 
to calculate the SC and HC are obtained by this simulated 
PCA spectrum. Aql X-1 SC and HC thus obtained are nor- 
malised by the average Crab SC and HC respectively. In 
figure [2] the point marked by an open plus sign with error 
bars represents this normalised Aql X-1 SC and HC. Thus 
we conclude that during the first Suzaku observation, Aql X- 
1 was in the BS. Henceforth this observation will be referred 
to as the BS observation. 

The PCA spectra from simultaneous observations (Ta- 
ble [l]) are used for comparison with the Suzaku spectra. We 
have used the Std2 spectra and response files provided in 
the standard products. Energy channels were appropriately 
regrouped before the final spectral fitting. 

2.2 Suzaku 

For each front-illuminated XIS detector (XISO and XIS3) we 
extracted the source spectra using a 260 arc-sec circular ex- 
traction region centered on the source. The energy scale is re- 
processed using the xispi task. Background spectra are ex- 
tracted using appropriate circular regions outside the source 
region. Corresponding response files are generated using the 
xisrmfgen and xissimarf gen tools. The spectra from re- 
spective XISO and XIS3 detectors are then added together 
using the addascaspec tool. 

During high count rates, as is the case for the BS ob- 
servation, there is a possibility of pile-up in the XIS detec- 
tors. To check for pile-up we extracted the source spectrum 
from an annular region with inner radius of 40 arc-sec and 
an outer radius of 260 arc-sec. Corresponding response files 
were also generated. Figure [3] shows the ratio of XISO back- 
ground subtracted spectrum extracted from the full circular 
region and from the annular region. The constant ratio in- 
dicates that both these spectra are similar and differ only 
in the observed count rates. Similar result is seen for XIS3 
detector also. This confirms that the BS observation is not 
affected by pile-up. 

The data processing version used for the HXD cleaned 
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Figure 2. Figure shows Colour-Colour diagram of Aql X-1 
for the 2007 outburst. The open circles with error bars are 
points during the peak of the outburst where no Suzaku ob- 
servations are available. Points marked with filled circles, filled 
stars and triangles are from RXTE observations simultaneous 
to the second, third and fourth Suzaku observations respec- 
tively. Note that there are no RXTE observations close to the 
first Suzaku observation. Points marked with rectangles are 
remaining RXTE observations after the decay of the outburst 
where no simultaneous Suzaku observations are available. The 
point marked with an open plus sign and error bars is obtained 
from the PCA spectra simulated using the best fit model for 
the first Suzaku observation. Error bars on points representing 
the EIS of Aql X-1 are suppressed for clarity except one point 
in the upper left hand corner. 
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Figure 3. Figure shows a plot of the ratio of XISO spectrum 
for BS observation, obtained using a full circular source region 
of radius 260 arc-seconds to the spectrum obtained from an 
annular region of inner radius 40 arc-seconds and outer radius 
of 260 arc-seconds at different energies. The constant ratio 
implies that both these spectra differ only in the observed 
count rates and there is no pile-up at the centre of the CCD 
during the BS observation. 



products is 2.1.6.15. Hence the PIN spectra are extracted 
using the cleaned event files following the standard anal- 
ysis threads on the Suzaku websiteQ The PIN non- X-ray 
background is extracted from the observation specific model 
provided by the instrument team. It is combined with the 
standard model for cosmic X-ray background. Appropriate 
response files for XIS nominal pointing available from the 



^ http:/ /heasarc. nasa.gov/docs/suzaku/analysis/abc/ 
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Table 2. XSPEC Model Components used 



Model 


XSPEC Components 




Ml (a/b) 
M2 


C^ X wabs (diskbb + diskline H 
Ct X wabs (BBody/BBodyrad 


- nthcomp)"'* 
■f nthcomp) 



f: Constant used to get relative normalisation between XIS 
and HXD spectra 

a: The input photons to nthcomp are from accretion disk 

b: The input photons to nthcomp are from an hidden 

blackbody emitting NS surface 

calibration database CALDB (version XIS-20101108, HXD- 
20101202) are usedQ 

The XIS spectra are grouped such that the final spectra 
have ~ 250 energy channels so that there are about three 
channels per energy resolution of ~ 150eV. Energy channels 
corresponding to the energy range of 1.75-1.95 keV are ig- 
nored since the response files do not adequately remove the 
absorption edge present in this energy range. Ten consecu- 
tive energy channels are grouped into one energy channel of 
the final grouped PIN spectrum. The net count rates for the 
XIS and PIN spectra are given in Table [U We use XSPEC 
V12.5.1n to fit the spectral models to the observed X-ray 
spectrum. 



3 RESULTS 

We have used two different phenomenological models to de- 
scribe the observed spectra of the source in both the EIS 
and BS. In the first model (henceforth referred to as Ml) 
we model the soft component using multi-temperature disk 
blackbody (DEB) and the hard component as Comptonised 
emission from the boundary layer. The hot boundary layer is 
assumed to completely cover the neutron star surface. This 
Eastern-like model has two variants. One, where the input 
photons to the boundary layer are from the accretion disk 
(Mia) and second where the input photons are from the 
underlying blackbody emission of the neutron star surface 
(Mlb). The presence of a standard accretion disk allows for 
the possibility of a broad Iron line emission in this model 
which we include using the XSPEC model "diskline". In the 
second Western-like model (M2) we model the soft com- 
ponent using a single-temperature blackbody (BE) emitted 
from the boundary layer. The hard component is modelled 
as Comptonised emission from an inner hot disk where the 
input photons are from the boundary layer. 

The Comptonised spectra for both the mod- 
els is represe nted by the XSPEC func tion "nth- 
comp" fiZdziarski. Johnson fc Magdziar j 1 19961 : 



IZvcki. Done fc Smithll 19991 ). We apply this Comptonization 
model mainly because it can have input seed photon 
populations of different nature and distribution (pure 
blackbody, multicolour disk) thus allowing to test different 
emission geometries. In the Mia model the input photon 
distribution is due to disk blackbody and hence the input 
temperature is set equal to the maximum temperature 

^ http:/ /heasarc. gsfc.nasa.gov/docs/heasarc/caldb/suzaku/docs/ 
suzaku_caldbhistory. html 



of the disk blackbody emission. In the M2b model, input 
photon distribution is due to the blackbody emission from 
the NS surface and hence the input temperature is a free 
parameter. On the other hand in the M2 model, the photon 
distribution is due to neutron star boundary layer and 
hence of a blackbody shape and the input tempertature is 
set equal to the temperature of this neutron star boundary 
layer. 

A multiplicative component ("wabs") is introduced in 
all spectral fitting to model the absorption of the source 
spectrum by intervening material along the line sight. We 
note here that the average Dickey & Lockman nH value in 
the direction of Aql X-1 is 0.33 x 10^^ cm"^. The nH val- 
ues obtained in our fits agree with this average value. An- 
other constant multiplicative term is used in both the models 
for the relative normalisation of the XIS and corresponding 
HXD spectra. The cross-calibration constant for XIS nomi- 
nal pointing as suggested by the Suzaku team is 1.16±0.014 
30. The value for the constant that we get from our fits agree 
with those suggested by the Suzaku team within errors. The 
XSPEC model components used for M2 and Ml are tabulated 
in Table O While fitting the near simultaneous RXTE spec- 
tra, we fix the absorbing column density, nH to the corre- 
sponding uH value derived from the Suzaku spectral fit since 
RXTE data alone cannot constrain the nH value. 

The normalisation parameter of the XSPEC BBody com- 
ponent gives the blackbody luminosity Lbb and the normal- 
isation parameter of the XSPEC BBodyrad component gives 
the blackbody radius Rbb- Similarly from the parameter 
values of the XSPEC nthcomp component we can estimate the 
total Compton luminosity Lcomp, the required luminosity of 
the input photons Li„p, the amplification factor (A) and 
when kTbb is a free parameter the radius of the blackbody 
Rbb- From the XSPEC diskbb model we calculate the inner 
accretion disk radius {Ri„) and the disk luminosity Ldbb- 
The calculated values of these parameters allow us to con- 
strain the two models considered here. Note that throughout 
this paper for calculating luminosities we use distance to the 
source as D = 5.25 ± 1.25 Kpc (Rutledge et al. 2001) and 
the error quoted for luminosities take into account this un- 
certainty in the distance. 

The results from the observations Ell, EI2 and EI3 
when the source was in the EIS will be presented together 
in the next sub-section and that of the BS in the subsequent 
one. 



3.1 Spectral analysis of Aql X-1 during the 
Extreme Island State 

We start with analysing the Suzaku spectra using Model 
Mia. If the Iron line emission is omitted the best fit x^/'^ 
obtained are 284.50/241, 295.74/237 and 305.31/241 for ob- 
servations Ell, EI2 and EI3 respectively. Adding a Gaussian 
line at 6.4 keV improves the fit considerably. The width of 
the Caussian turns out to be > 0.5 keV, indicating the line is 
broad. Hence, we adopt the more physical "diskline" model 



http: / /www. astro. isas.jaxa.jp/suzaku/doc/suzakumemo/suzaku 
memo-2007-ll.pdf 

http:/ /www. astro. isas.jaxa.jp/suzaku/doc/suzakumemo/suzaku 
momo-2008-06.pdf 
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Table 3. Parameter values with respective la errors for model Mia during the EIS 



Mia Parameter 




Ell 




EI2 




EI3 




402053020 


93405-01-04-01 


402053030 


93405-01-05-01 


402053040 


93405-01-06-00 


nH(1022)cm-2 


0.311 ±0.008 


0.311{frozen) 


0.287 ±0.007 


0.287(frozcn) 


0.445 ± 0.009 


0.445(frozen) 


T,„{keV) 


0.85 ±0.02 


0.85 ±0.11 


0.82 ±0.02 


0.97 ±0.09 


0.82 ±0.02 


0.79 ±0.05 


iDSs(1035erg/s) 


5.1 ±2.4 


5.0 ±2.5 


4.7 ±2.2 


6.9 ±3.2 


5.2 ± 2.5 


6.9 ±3.3 


«L(km) 


6.6 ±0.2 


6.8 ± 1.1 


6.7 ±0.2 


5.9 ±0.9 


7.1 ±0.2 


9.0 ± 1.1 


7 


1.71 ±0.04 


1.85 ±0.05 


1.78 ±0.03 


1.73 ±0.04 


1.71 ±0.05 


1.74 ±0.04 


kTe(keV) 


20.0{frozen) 


20.0(frozen) 


20.0(frozen) 


20.0(frozen) 


13.6t^;^ 


9 7+2-9 
'-1.7 


icompClO^^erg/s) 


25.5 ± 12.1 


25.4 ± 12.1 


29.4 ± 14.0 


30.6 ±14.5 


21.8 ± 10.4 


19.8 ±9.4 


iinp{1035erg/s) 


4.6 ±2.2 


6.1 ± 2.9 


6.1 ± 2.9 


6.1 ± 2.9 


4.3 ±2.0 


4.6 ±2.2 


A 


5.54 


4.16 


4.82 


5.01 


5.07 


4.30 


CX 


1.05 ±0.06 




1.11 ±0.05 




1.10 ±0.07 






266.70/240 


30.48/28 


260.93/236 


37.03/28 


287.22/240 


24.23/23 



t Propagating distance errors in calculating error on i?i„ would increase the errors to ~ 5 km for all observations. 
I C is the cross-calibration constant between the XIS and HXD spectra. 
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Figure 4. The left hand panel shows the Suzaku and the right hand panel shows the corresponding RXTE-PCA unfolded spectrum 
for the EI3 observation. Observed data points are marked in black. In the left panel points below 10 keV are from the XIS detector 
and the points marked with a circle are from the HXD detector. Red line shows the contribution from the disk blackbody component, 
the Blue line shows the contribution of the Comptonised component and the green line shows the contribution from the Fe line. 
The bottom panels show the residues as a function of energy with (middle panel) and without (bottom panel) the Fe line. 



to represent the Iron line. However, the data is not good 
enough to constrain all the parameters of the model and 
hence we fix some of them (emissivity index /3 = —2.0, in- 
clination angle i = 45°, inner radius Rin = lOGM/c^ and 
outer radius Rout = lOOOGM/c'^). Thus only the normali- 
sation of the line is a fitted parameter. The same model is 
applied to near simultaneous RXTE observations. Removing 
the Iron line component for the RXTE spectra, gives best 
fit x^/^ values of 38.39/29, 70.08/29, 64.21/24 for Obs. Id. 
93405-01-04-01/05-01/06-00 respectively, which are signifi- 
cantly higher than obtained with the line feature. Thus, 
the Iron line feature is detected both in the Suzaku as well 
as in the RXTE spectra. 

The electron temperature of the Comptonizing bound- 
ary layer component, kTin is unconstrained for observations 
Ell and EI2 for both Suzaku and RXTE data with a lower 
limit of above 15keV. Hence we fix kTin = 20 keV for these 
observations. The spectral parameters and the best fit /v 
obtained for the Suzaku and RXTE spectra are given in Ta- 
ble [S] Figure |4] shows the E^F{E) spectrum for Suzaku ob- 



servation EI3 and corresponding RXTE spectrum along with 
the Mia model components. Note that the best fit spectral 
parameters are consistent within errors for the Suzaku and 
RXTE observations for all three sets. This similarity and 
consistency between the two different instruments, shows 
that the spectral extraction and response generation is ro- 
bust and correct. 

The inner radius of the accretion disk turns out to be 
~ 7 km which is rather small since it needs to be larger than 
the radius of the neutron star Rs 10 km. However, this 
may be due to the uncertainty in the distance to the source. 
Assuming that the distance is 6.5 instead of 5.25 kpc used 
here would lead to physically acceptable inner radii of > 10 
km. Moreover, the colour factor assumed here is f — 1.7 and 
a slightly higher value could also correct this discrepancy. 
Note that the luminosity of the disk is of the same order as 
the input luminosity into the Comptonizing medium. This 
is an important consistency check since the model assumes 
that input photons arise from the accretion disk. 

Next we consider the possibility that the input photons 
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to the Comptonizing medium arise from a blackbody emis- 
sion of the underlying neutron star surface (Model Mlb). 
Thus in this model temperature of the blackbody emission, 
kThb is an extra free parameter. The best fit /v values are 
268.35/239, 266.18/235 and 260.83/239 for Ell, EI2 and EI3 
observations respectively. The RXTE spectral fits are also 
consistent with the Suzaku fits, although the extra freedom 
of kTtb does not allow the spectral parameters to be well 
constrained. But it is noted that for Ell and EI2, the val- 
ues of kTbb and Ti„ coincide within errors whereas for EI3 
kTbb < Tin which is not compatible with the geometry of this 
model. The calculated inner disk radii values for this model 
turn out to be larger at ~ 9 km. However, the radius of 
the blackbody assumed to be from the neutron star surface, 
turns out to be rather small (~ 5 km) for the first two ob- 
servations, namely Ell and EI2. Again, physical consistency 
may be obtained by using a larger distance to the source 
and/or including a colour factor in the blackbody radius cal- 
culation. Also for the EI3 observation, the obtained radius 
of the blackbody is Rbb ~ 18±5 km which is greater than 
the calculated inner disk radius Ri„ = 7.5 ± 0.3 km. There- 
fore, perhaps the more physical picture for EIS is that the 
input photons for Comptonization are mainly contributed 
by the accretion disk making Mia model more feasible than 
the Mlb model in this case. 

We consider next the alternate geometry where the 
Comptonised component arises from an inner hot disk and 
the boundary layer produces a blackbody emission (Model 
M2). When this model is used to fit the three EIS spec- 
tra the best fit jv for observations Ell, EI2 and EI3 
are 300.87/241, 306.83/237 and 267.50/241 respectively. In 
this model the input photon temperature for Comptoniza- 
tion component kTbb is kept equal to the corresponding 
blackbody temperature kT . The calculated blackbody ra- 
dius Rbb for the three observations are 3.6 ± 0.9, 3.5 ± 0.8 
and 4.1 ± 1.0 km respectively. These values are significantly 
smaller than the neutron star radius of ~ 10 km. The values 
can be made consistent only by assuming a large distance to 
the source of > 11 kpc or by assuming an unphysical large 
colour factor of ~ 3. Note that for a radius of 3 km the 
emitting region is only one-tenth the neutron star surface. 
This makes the geometry unfavourable as compared to the 
Ml model. 

Moreover, the values obtained are significantly larger 
than the ones obtained for Ml models. Even if we include 
a broad Iron line (which in the absence of a cold disk may 
not be physical) the x^/j^ values are 293.99/239, 285.49/235 
and 260.46/239 which are higher than the ones obtained 
for the Ml model. Consistent results are obtained from the 
RXTE data which provide jv of 34.40/28, 87.43/29 and 
24.33/23 respectively, again significantly higher than what 
was obtained earlier. 

3.2 Spectral analysis of Aql X-1 during the 
Banana State 

We now come to the spectral analysis of the first Suzaku 
observation taken during the decline of the 2007 Aql X-1 
outburst. During this observation the overall X-ray source 
flux is higher than that during the next three observations 
when the source is in the EIS. It was seen that for this BS 
Suzaku spectrum there is an excess of count rate at the high 
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Figure 5. The Mia model spectral components (listed in Ta- 
ble[4]along with the unfolded spectrum for the BS Suzaku ob- 
servation. Black points below 10 keV show the XIS spectrum 
and black points marked with circles show the HXD spectrum. 
The multicoloured disk blackbody component is shown in red, 
the Comptonised component is shown in blue, the power-law 
component in magenta and the Fe line component in green. 
The bottom panels show residues as a function of energy with 
(middle panel) and without (bottom panel) Fe line. 



energy end of the spectrum which could not be modelled 
by thermal Comptonization. We have added a hard power- 
law component (with photon index F = 0) to all spectral 
fitting to account for the high energy excess. We note here 
that energy contribution of this component diverges with 
energy and hence there should be a cutoff at both the high 
and low energy end. But for the given data neither the pho- 
ton index nor the cutoffs at low and high energies could be 
constrained. Therefore luminosity contributed by this com- 
ponent has been calculated only in the energy range of 0.7- 
70.0 keV in which the spectral fitting has been performed. 

Table|4]lists the best fit parameters for the variations of 
the cold accretion disk model (Mla/b) while Figure [5] shows 
the unfolded spectrum with the different components for the 
Mia model. The electron temperature of the Comptonizing 
region is smaller (fcTe ~ 2.3 keV) as expected and the in- 
ner radius of the accretion disk is significantly larger at ~ 40 
km. As in the EIS fitting, the disk luminosity is nearly equal 
to the input luminosity (as expected in this geometry), al- 
though the luminosities have increased by nearly an order of 
magnitude. Here, in Mlb model, the radius of the underlying 
neutron star, Rbb ~ 10 km, which is larger (and hence more 
physical) than that obtained for the EIS state. Note that the 
high energy power-law component contributes around 1.5% 
of the total luminosity. Without this component the /'^ is 
282.81/234 for Mia model and 279.90/233 for M2b model. 
For both Mia and Mlb model, the emissivity index of the 
"diskline" model can be constrained to g ~ —3. 

The alternate geometry of an hot inner disk (M2) pro- 
vides an unacceptably larger x^/'^ of 359.99/235. The inclu- 
sion of a broad Gaussian for the Iron line (which may not 
be physical given the absence of a cold disk) reduces it to 
338.43/233, which is still significantly higher than what is 
obtained for the Ml models. However, in this case the radius 
of the blackbody emitting boundary layer turns out to be 
~ 18 km and hence is physically possible. 
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Table 4. Parameter values of model Ml(a/b) for Obs Id 
402053010 



the optical depth of the medium, for a spherical geometry 
can be estimated by l|Sunvaev fc Titarchuklll980l 'l 



IVL _L Jr^circLIIlStGr 


iVi ±ci V di Ilea 


A/TI K \/ciliioc; 
iVl ± U V di Llco 


nH(1022) cm-2 


0.370 ± 0.009 


0.366 ± 0.008 


T™(keV) 


0.66 ± 0.02 


0.71 ±0.02 


Ldbb (erg/s) 


c; o 1 o Q\ \y 1 n36 
{o.o ± Z.o) X iU 


f '] n Q J- A n\ \y 1 n36 
t^iU.o ± 4.yj X iU 


R.n (Km)" 


37.1 ± 1.3 


43.1 ± 1.7 


pi, 


O.O(frozen) 


O.O(frozen) 


LpL (erg/s) 


(2.6 ± 1.2) X 10^5 


(2.7 ±1.3) X 10^^ 


7 


1.94 ± 0.09 


9 1+0.6 
^•^-0.3 


kTe (keV) 


2.3 ±0.1 


2 4+°-2 


kTtt (keV) 




1 1+0.3 
^•^-0.4 


Lcomp (erg/s) 


(10.6 ± 5.0) X 1036 


(6.1 ±2.9) X 10^6 


Linp (erg/s) 


(4.6 ±2.2) X lO^^* 


(4.3 ±2.1) X lO^'^ 


A 


2.30 


1.41 


Bbb (Km)'* 




9.5 ±6.9 


q 


o 1+0.6 
■^■^-1.7 


_o 1+0.7 
'3-i_4.7 


c 


1.19 ±0.09 


1.16 ±0.08 




246.78/233 


246.10/232 



a: Propagating error on distance increases error in Ri„ to 25 km 
for Mia and to 29 km for Mlb. 

b: The power law should have both a low and high energy cut-off 

which cannot be constrained by the present data. 

c: Luminosity due to the powerlaw tail is calculated in the energy 

range 0.7-70.0 keV. See results section for more details. 

d: Propagating error on distance increases error in Rbb to 7.2 km. 



4 SUMMARY AND DISCUSSION 

In this work, we analyse Suzaku and RXTE observations of 
Aql X-1, to constrain the geometry of the source in both the 
low flux EIS and high flux BS, and to report the discovery 
of a hard X-ray tail in the BS. In particular, a standard ac- 
cretion disk with a hot Comptonizing boundary layer whose 
input photons are from the disk is the preferred geometry for 
both states. The alternate model of a hot Comptonizing in- 
ner disk and a blackbody emitting boundary layer gives sig- 
nificantly larger for all four Suzaku observations includ- 
ing three nearly simultaneous RXTE observations. More- 
over, for the EIS observations, the radius of the blackbody 
emitting boundary layer turns out to be unphysically small, 
~ 3 km. The presence of a broad Iron line feature provides 
further support for the existence of a cold accretion disk in 
the source. The preferred geometry is the same as i nferred 
from broad band study of many other NS-LMXBs (|Barretl 
I2OOII ) and hence seems to be general. 

For all observations, the luminosity of the input pho- 
tons to the Comptonizing hot boundary layer is similar to 
that of the disk and hence it is possible that the main source 
of input photons is the disk itself. While it may also be pos- 
sible that the source of soft photons could be the underlying 
blackbody emitting neutron star, statistically the data can- 
not differentiate between these scenarios. 

Our results confirm that the electron temperature of the 
Comptonizing boundary layer in EIS is higher as compared 
to same in the BS. Also the inner radius of the accretion disk 
is higher during the BS observation. The spectral slope 7 ^ 
1.8 remains nearly constant, implying that the optical depth 
has increased in the high flux state. In the "nthcomp" model. 



^ thcomp — 



9 1 

4 ^ (fcTe/mec2)r(l + r/3) 



1/2 



(1) 



The optical depth r ~ 2 for the observations Ell and EI2 
and ~ 5 for the EIS observation. In the BS r ~ 13. The 
difference in the electron temperature is the primary cause 
for the two states to have different locations in the colour- 
colour diagrams. A larger number of observations, sampling 
the source at different times of the outburst will provide 
more information regarding the complete temporal evolution 
of the source. 

For the EISs, when the input photons are assumed to 
be from the disk (Table [3]), the luminosity of the disk is 
~ 5 X 10''^ ergs s~^. The luminosity of the boundary layer, 
which is the luminosity in the Comptonised component mi- 



nus the input luminosity Lc 



2.5 X lO-"" ergs 



s~ . Thus energy released in the boundary layer is about 5 
times larger than that of the disk. When the input photons 
are taken to be the underlying neutron star surface, bound- 



ary layer luminosity, Lc 



2.8 X 10"^^ ergs s"^ is about a 



factor of four more than the disk luminosity ~ 7 x 10 ergs 
s~^. For the BS this ratio is ~ 1 when the input photons are 
disk photons and ^0.7 when they arise from the underlying 
star surface. Theoretically, the ratio of the energy output in 
the boundary layer to that of the disk can range from 1-3 
dependin g on the spin and equatiori of state of the neutron 
star (e.g. iBhattacharvva et aUl2000l ). A ratio of five would 
imply an extremely soft equation of state. However, the sim- 
ple phenomenological models used here do not warrant one 
to make concrete statements but one can conclude that the 
boundary layer is more luminous than the disk, which is 
consistent with the expected energy release ratio. 

An interesting result is the detection of an hard 
X-ray power-law emission in the high flux state which 
is only 1.5% of the total source luminsoity. Such hard 
tails have been detected in Z sou rces and only one Atoll 
source GX 13+1 earlier jPaizis et al. 2006,). The origin of 
this emission is unknown a nd maybe due to some non- 
thermal emission in a jet ()Migliari fc FendeiJ 120061 ). In- 
deed, a weak radio emiss ion has been detected in the high 
soft state of Aql X-1 (|Rupen. Mioduszewski. fc Dhawar] 
2004). which is also o bserved in other atoll sources 
(jMigUari fc Fender] |2006| ). Alternatively, the hard X-ray 
emiss ion could be due to an hybrid non-thermal emission 
(e.g. [Coppil Il999|) or bulk motion Co mptonization (e.g. 
iTitarchuk. Mastichiadis fc Kvlafis|[l997l ). Multiple observa- 
tions of the source may reveal whether the hard X-ray tail is 
variable and whether it is correlated with the radio emission. 

The work highlights the utility of broad band spectral 
analysis of LMXB. Apart from other Suzaku observations 
the upcoming Indian satellite ASTROSAT will be useful in 
such studies of LMXBs containing weakly magnetised neu- 
tron stars as it will provide good timing and spectral reso- 
lution. 
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